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The method of preliminary group classification is rigorously defined, enhanced and related 
. to the theory of group classification of differential equations. Typical weaknesses in pa- 

pers on this method are discussed and strategies to overcome them are presented. The 

■ preliminary group classification of the class of generalized diffusion equations of the form 
^ I ut = f{x,u)u'^ + g{x,u)uxx is carried out. This includes a justification for applying this 

method to the given class, the simultaneous computation of the equivalence algebra and 
equivalence (pseudo)group, as well as the classification of inequivalent appropriate subal- 
gebras of the whole infinite-dimensional equivalence algebra. The extensions of the kernel 
algebra, which are induced by such subalgebras, are exhaustively described. These results 
i-S^ . improve those recently published in Commun. Nonlinear Sci. Nunier. Simul. 

^ ; 1 Introduction 

Group classification of differential equations is an efficient tool for investigating symmetry prop- 
erties of classes of differential equations. These are differential equations that include arbitrary 
^ ■ constants or functions of the independent and dependent variables as well as of derivatives of 

! the dependent variables up to a certain order. It is known for a long time that depending on 

the value of these arbitrary elements the resulting differential equations from the given class can 

■ have different Lie invariance groups. The first examples of group classification were presented by 
. Sophus Lie for the class of second order linear partial differential equations [15] and the class of 

second order ordinary differential equations [16]. Later, Ovsiannikov began the rigorous devel- 
opment of the theory of group classification [28]. In short, the solution of the group classification 
problem consists in finding the kernel of Lie invariance groups (i.e. those Lie symmetries that 
are admitted for all values of the arbitrary elements) and all inequivalent extensions of Lie in- 
^ ' variance groups with respect to the kernel group. The equivalence involved means the similarity 

. of equations up to transformations from a certain equivalence group (e.g. usual, generalized or 

conditional equivalence), see [32] for more detailed information. 

For classes of differential equations being of simple structure (e.g., ones parameterized only 
by constants or functions of the same single argument), the corresponding group classification 
problems can be completely solved via compatibility analysis and explicit integration of the 
determining equations for Lie symmetries depending on values of the arbitrary elements and 
up to the equivalence chosen. Complete group classification can also be carried out for classes 
of differential equations possessing the normalization property. The algebraic method of clas- 
sification effectively works for such classes. See the next section and also [32, 37] for a more 
comprehensive review on different methods of group classifications. 

In the situation where the class depends in a more complicated way on its arbitrary elements, 
it may happen that both the determining equations are too difficult to be directly solved and the 
application of the algebraic method does not give the exhaustive solution. In this case, however, 
at least a partial solution of the group classification problem, known as preliminary group clas- 
sification, is possible. The basic idea of preliminary group classification is to study only those 
extensions of the kernel group that are induced by the transformations from the corresponding 
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equivalence group. The problem of finding inequivalent cases of such Lie symmetry extensions 
then reduces to the classification of inequivalent subgroups (resp. algebras) of the equivalence 
group (resp. algebra). This approach was first described in [1] and became prominent due to 
the paper [11]. 

Despite the approach of preliminary group classification is rather common, it is not well de- 
veloped up to now. The basic mechanisms were formulated in [11] as two propositions for the 
specific class of nonlinear wave equations of the form vu = f{x, Vx)vxx + g{x, Vx) and were later 
adopted in other papers for respective classes of equations. In the present paper, we state a 
stronger version of these propositions for general classes of differential equations. Another weak- 
ness commonly observed is that, when the equivalence algebra 0~ of the class of equations under 
consideration is infinite dimensional, only Lie symmetry extensions induced by subalgebras of a 
finite-dimensional subalgebra of 0~ are investigated, without giving any sound justification 
for the choice of Qq . In fact, this restriction is needless as it is possible to classify subalgebras of 
infinite-dimensional algebras in much the same way as subalgebras of finite-dimensional algebras 
[3, 4, 5, ?, 9, 14, 18, 31, 32, 38]. It can even be simpler to classify low-dimensional subalgebras 
of the whole infinite-dimensional equivalence algebra q"" as the adjoint action related to g~ is 
more powerful and allows for greater simplification than the adjoint action corresponding to 
the finite-dimensional subalgebra . One more common weakness in papers on the subject is 
that usually only extensions induced by one-dimensional subalgebras of equivalence algebras are 
studied. Moreover, these one-dimensional subalgebras (of a finite-dimensional subalgebra 
of 0~) are classified only with respect to the restricted equivalence relation which is generated 
by the adjoint representation of . This leads to an overly large number of inequivalent sub- 
algebras compared to the list of one-dimensional subalgebras that would be obtainable if the 
classification was done using the adjoint representation of the entire 0~. 

In the present paper, we comprehensively carry out preliminary group classification for the 
class of (1 -|- l)-dimensional second order quasilinear evolution equations of the general form 



where / and g are arbitrary smooth functions of x and u, and g ^ 0. The class (1) was considered 
in the recent paper [20] but results obtained therein are not correct. It is reviewed above that 
there are a number of typical weaknesses in papers on preliminary group classification, and 
results of [20] properly illustrate these weaknesses, cf. the first paragraphs of Sections 4 and 6 and 
Remark 8. This is why we aim to accurately present the revised preliminary group classification 
of the class (1) and to give all calculations in considerable detail. 

The class (1) was considered in [20] as a class of generalized Burgers equations as it includes 
the potential Burgers equation as a particular element for the choice f = g = 1. This class also 
contains (1 + l)-dimensional linear evolution equations, which correspond to the values / = 
and g not depending on u. As a prominent example for a linear differential equation, one can 
recover the linear heat equation by choosing / = and g = 1. An important subclass of the 
class (1) is the class of (1 -I- l)-dimensional nonlinear diffusion equations of the general form 
ut = {F(u)ux)xi where F ^ 0. It is singled from the class (1) by the constraints gx = 
and / = gu- Moreover, any equation of the form (1) with fx = gx = Ois reduced to a diffusion 
equation by a simple point transformation acting only on the dependent variable u. The solution 
of the group classification problem for this class by Ovsiannikov [27] (see also [1, 28]) gave rise 
to the development of modern group analysis. 

The class (1) is included in the wider class of equations ut = F{t, x, u, Ux)uxx + G{t, x, u, Ux), 
for which the complete group classification was carried out in [3] by a method similar to that 
applied in the present paper. The fact that this class is normalized (cf. Section 2) played a 
crucial role in the entire consideration in [3]. However, as this class is essentially wider than the 
class (1), the corresponding equivalence algebras are rather different. This is why the results 
of [3] cannot be directly used for deriving the group classification of the class (1). 
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It should also be stressed that it is not natural to exclude linear differential equations from the 
present consideration. In fact, there are equations in the class (1) which are linearized by point 
transformations from the equivalence group C of this class. The most prominent example of 
such a transformation in the above class is the transformation of the potential Burgers equation 
to the linear heat equation by means of the point transformation u = [25, p. 122]. That 
is, u is a solution of the potential Burgers equation whenever « is a solution of the linear heat 
equation. In the course of preliminary group classification of the class (1) we encounter other 
examples of linearizable equations. Furthermore, the equivalence algebra of the subclass 
of (1), which is compliment to the subclass of linear equations and, therefore, singled out by the 
constraint + Qu 7^ 0, is much narrower than the equivalence algebra 0~ of the entire class (1). 
More precisely, the algebra gQ is singled out as a subalgebra of 0~ by the constraint huu = 0, 
cf. Theorem 1. 

The further organization of this paper is the following. The subsequent Section 2 discusses 
the theory of preliminary group classification. We generalize and extend assertions presented 
in [11] and formulate them rigorously using the modern language of group analysis. In Section 3 
we derive the determining equations for Lie point symmetries of equations from the class (1) and 
find the corresponding kernel of Lie invariance algebras. The equivalence algebra and the 
equivalence group C of the class (1) is computed in Sections 4 and 5, respectively. Throughout 
the paper, by the equivalence group we mean the Lie pseudo-group of point equivalence trans- 
formations (i.e., local equivalence diffeomorphisms), cf. [26] and references therein for theory of 
pseudo-groups. The reason for carrying out preliminary group classification is elucidated. In 
Section 6, we classify inequivalent one- and two-dimensional subalgebras of the essential subalge- 
bra of . The corresponding inequivalent cases of symmetries extensions of the kernel algebra 
are presented in Section 7 and supplemented with three- and four-dimensional extensions via the 
classification of all appropriate subalgebras of g"". The paper concludes with a short summary 
and further comments in Section 8. 



2 Enhanced method of preUminary group classification 

By now, the method of preliminary group classification was neither explained for general classes 
of differential equations nor properly related to the general group classification problem. This 
should be done first in this section before we study the preliminary group classification of (1). 
For this aim, we need a few notions of the theory of group classifications, which can be found 
in the recent paper [32]. 

The most essential notion concerns the formal definition of classes of differential equations. 
In general, a class (of systems) of differential equations is given by a system of / differential 
equations of the form L{x, U(^p^, 9{x, U(p))) = in m dependent variables u = (u^, . . . , n"*) and n 
independent variables x = (xi, . . . , Xn), where M(p) denotes the set of u's and all their derivatives 
up to order p. The differential functions 9{x, 'U(p)) = {0^{x, u^p-j), . . . , 6^{x, ii(p))) denote a tuple 
of k arbitrary elements that parameterize the given class of differential equations. The tuple 9 
is usually constrained to satisfy a system S of auxiliary conditions, 5(x, U(p), ^(g)(x, M(p))) = 0, 
in which x and ii(p) are regarded as independent variables. The set of solutions of this auxiliary 
system will also be denoted by S. In addition, this set can be further constrained by satisfying 
one or more nonvanishing conditions S(x, 0(q)(x, / 0. Putting together all these 

notions, we denote the class of differential equations with the arbitrary element running through 
the set S by C\s. The single elements of this class are denoted by Cq, respectively. 

Specifically, for the class (1) we have 6 = {f,g), and the arbitrary elements / and g depend 
only on x and u. Therefore, the associated auxiliary system S is formed by the equations 

ft fut fux futt futx fuxx 0, 

9t = 9ut = 9ux = 9utt = 9utx = 9uxx = 0- 
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The auxiliary conditions ft = and gt = play a special role. All the other auxiliary conditions 
can be taken into account implicitly. The nonvanishing condition associated with the class (1) 
is (7 7^ 0, i.e., we have = g. The condition g ^ should be explicitly included in the definition 
of the class (1) since equations of the same form with g = are of another (first) order, possess 
completely different transformational properties and are not related to equations with g ^ hy 
point or other reasonable transformations. 

Having properly defined classes of differential equations, it remains to introduce the notion 
of admissible transformations and normalized classes of differential equations in order to explain 
the general strategy of (preliminary) group classification. 

Definition 1. The set of admissible transformations in the class C\s is given by T{C\s) = 
{{6,9,ip) \ 8,9 £ S,{p T{6,6)}, where T{6,9) denotes the set of point transformations that 
map the system Cg to the system Cg. 

The set of admissible transformations can be used for many issues related to the problem of 
group classification. It can be considered as an extension or generalization of the equivalence 
group of a class of differential equations. Indeed, the usual equivalence group C of a class 
C\s is naturally embedded in the set of admissible transformations. In particular, it is given 
by the admissible transformations {9, ^9, ^\(x,u))j where $ is an equivalence transformation, i.e. 
\/9 & S: ^9 £ S. In this last tuple, ^\(x,u) £ T(^,<I>^) denotes the projection of $ to the space 
of variables x,u. The maximal point symmetry group Gq of the system Cg coincides with the 
set of admissible transformations from Cg to itself, i.e., Gg = T{9,9). 

Important properties of classes of differential equations, relevant for the problem of group 
classification, are given by different kinds of normalization with respect to point (resp. contact) 
transformations [30, 32]. 

Definition 2. The class of differential equations C\s is normalized in the usual sense if any 
admissible transformation is induced by a transformation of the (usual) equivalence group, i.e. 
y{9,9,^), 3^eG-: 9 = ^9 and yD = 

Denote by Qg the maximal Lie invariance algebra of the equation Cg. Using the above 
notations it is possible to obtain the general picture of the group classification problem. The 
first step in order to carry out group classification is the determination of the kernel (i.e., 
intersection) of maximal Lie invariance algebras of systems from the class £1^. The kernel is 
found by deriving the determining equations of Lie symmetries and splitting with respect to both 
derivatives with respect to u and the arbitrary elements 9. This gives those part of the maximal 
Lie invariance algebra Qg that is admitted for any value of 9. The subsequent step consists of 
determining the equivalence group G~ (resp. the equivalence algebra q^) of the class C\s. The 
equivalence group G~ is needed since it generates a natural equivalence relation on cases of 
symmetry extension of the kernel and hence they should be studied up to this equivalence. The 
final task is to describe all inequivalent cases of symmetry extension, i.e., values of 9 for which 

For the implementation of the above classification program, several special techniques have 
been developed. They either lead to the complete group classification or to a preliminary group 
classification of the given class. 

Complete group classification is often possible for normalized classes of differential equations. 
For such classes, symmetry extensions of the kernel algebra can only be induced by vector fields 
from the corresponding equivalence algebra. This reduces the group classification problem to 
the algebraic problem of classifying inequivalent subalgebras of the equivalence algebra. This is 
why we refer to this method as the algebraic method. Results on complete group classification 
of various classes of differential equations can be found, e.g. in [1, 3, 14, 31, 32, 33, 34, 38]. The 
equations studied in these papers all possess the normalization property. 
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Another method leading to the complete solution of the group classification problem consists 
of a compatibility analysis and direct integration of the determining equations of Lie symmetries 
of the given class [1, 2, 12, 21, 22, 23, 24, 28, 37]. It was indicated in the introduction that it is 
often only for rather simple classes that this method works. 

Complete preliminary group classification employs essentially the same techniques that are 
used for complete group classification within the framework of the algebraic method. The 
main difference is that the underlying class does not possess the normalization property. This 
implies the existence of extensions of the kernel algebra that are not induced by subalgebras of 
the equivalence algebra. In turn, for normalized classes of differential equations the results of 
complete preliminary group classification and complete group classification coincide [30, 32]. 

In most papers on preliminary group classification only a partial solution of the corresponding 
problems is achieved since usually not the whole equivalence algebra is used for an investigation 
of cases of symmetry extensions. This is why we refer to this method as the method of partial 
preliminary group classification. It is the most incomplete and heuristic method of group clas- 
sification, as there are often no obvious criteria which subalgebras of the equivalence algebra to 
single out for an investigations of symmetry extensions of the kernel algebra. Results on partial 
preliminary group classification are presented, e.g., in [1, 11, 20, 36]. 

On the side of complete group classification, the theoretical background was already set- 
tled [17, 28] and extended [30, 32]. It remains to detail the framework of preliminary group 
classification. In its essence, it rests on the following two propositions, which were first for- 
mulated without proof in [11] for the class of equations investigated. We present an enhanced 
version of these propositions for general classes of differential equations. 

Proposition 1. Let a be a subalgebra of the equivalence algebra 0~ of the class C\s, a C and 
let 6^{x,U(^j.^) £ S be a value of the tuple of arbitrary elements 6 for which the algebraic equation 
9 = 6^{x,U(^j.)) is invariant with respect to a. Then the differential equation Cgo is invariant with 
respect to the projection of a to the space of variables {x,u). 

Proof. Choose an arbitrary operator Q from a and consider the one-parameter group Gi gen- 
erated by this operator. As the equation 6 = 6^{x,U(^r^) is invariant with respect to Gi, any 
transformation T from Gi maps the corresponding equation Cgo from the class C\s to itself. 
This means that the projection PT of T to the space of variables {x, u) is a point symmetry 
of Cgo . Therefore, the projection PGi of Gi is a point symmetry group of Cgo and its generator, 
which is the projection of the operator Q, belongs to the Lie invariance algebra of Cgo. □ 

Proposition 2. Let Si be the subset of S that consists of all arbitrary elements for which 
the corresponding algebraic equations are invariant with respect to the same subalgebra of the 
equivalence algebra g'" and let be the maximal subalgebra of g~ for which Si satisfies this 
property, i = 1,2. Then the subalgebras ai and 02 are equivalent with respect to the adjoint action 
of G^ if and only if the subsets Si and S2 are mapped to each other by transformations from G^ . 

Proof. Assume that a2 = Xiai, where T G G~ and Tt, denotes the associated push-forward of 
vector fields. For 6^ G Si the algebraic equation ^ = is invariant with respect to ai. Since 
T is an equivalence transformation, we also have that TO^ G S. By supposition, the equation 
= is invariant with respect to T^oi = 02. This implies that TO^ G ^2 from which it can be 
concluded that TSi C ^2. Similarly, for 6^ £ S2, the algebraic equation 6 = 6^ is invariant with 
respect to 02 and T~^9^ G S. As T~^a2 = ai, the algebraic equation 6 = T~^6^ is invariant with 
respect to ai, which implies that T~^0^ G 5*1. From this last condition we obtain TSi D S2. 
It therefore can be concluded that there exists a bijection between Si and ^2, generated by a 
transformation from 

Conversely, suppose that ^2 = TSi for T G G^ . If = 6^ is invariant with respect to Oi then 
9 = T6^ is invariant with respect to T^ai. As 6^ is arbitrary, this implies that T^ai C a2. In a 
similar manner as in the previous paragraph, we can show that 7^~^02 C ai using the inverse 
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transformation of T. Then we have 7^T^ ^02 C T^ai and thus 02 C T^Oi. This is why 02 = T^ai 
must hold, which completes the proof of the proposition. □ 



Roughly speaking, the first proposition defines the method of how to construct cases of sym- 
metry extensions if the equivalence algebra of the class of differential equations to be investigated 
is already known. The second proposition then states that the problem of finding inequivalent 
cases of such symmetry extensions of the kernel algebra is reduced to the algebraic problem of 
the classification of subalgebras of the equivalence algebra. 

Remark 1. Within the set Si defined in Proposition 2, there is an equivalence relation generated 
by transformations from G~ whose push-forwards to vector fields preserve the subalgebra 
of Such transformations form the normalizer of the subgroup of G~ associated with aj. This 
equivalence relation can be used to choose simple forms of representatives of the set Si. 

This now completes the picture of the methods available for general group classification prob- 
lems. It should be clear that these methods apply to different classes of differential equations. 
This is why it is essential to investigate properties of the given class before choosing a particular 
method of group classification. This is done in the present paper. It is shown in the subse- 
quent sections that the class (1) is not normalized. Moreover, a compatibility analysis of the 
determining equations of Lie symmetries of this class is also an overly complicated task. This is 
why it cannot be expected to solve the complete group classification problem for (1) in a rea- 
sonable way. Still, the given class is adequate to be investigated using the method of complete 
preliminary group classification. 

Recall that, as mentioned in the introduction, the class (1) is contained in the wider class of 
equations of the general form ut = F(t, x, u, Ux)uxx + G{t^ x, u, Ux), which is normalized and for 
which the group classification problem was solved in [3]. 

The following folklore assertion is true. 

Proposition 3. The kernel (common part) = flegs ^/ maximal point symmetry 
groups Ge, 9 € S, of systems from the class C\s is naturally embedded into the (usual) equivalence 
group of this class via trivial (identical) prolongation of the kernel transformations to the 
arbitrary elements. The associated subgroup G^ of G~ is normal. 

Proof. Let To be an arbitrary element of G^, i.e. To is a point symmetry transformation for 
any equation from the class C\s. Denote by To the trivial prolongation of To to the arbitrary 
elements 0, To^ = 9. The transformation To obviously belongs to since it maps any equation 
from C\s the same equation in the new variables and therefore saves the entire class C\s. 

Taking an arbitrary transformation T E C, consider the composition T~^'TqT. In order 
to check that G^ is a normal subgroup of we should prove that this composition belongs 



to G^. We fix any 61 G 5 and denote T9 by 0. Then %T0 = 9 and hence T~^%T9 = 9. 



of G^. Therefore, T ^7ol~, which is the trivial prolongation of PT" ^ToT to the arbitrary 



Properties of G'~^ described in Proposition 3 were first noted in different works by Ovsiannikov 
(see, e.g., [29] and [28, Section II.6.5]). Another formulation of this proposition is given in [17, 
p. 52], Proposition 3.3.9. 

Corollary 1. The trivial prolongation of the kernel algebra to the arbitrary elements is 
an ideal in the equivalence algebra g~. 

By definition, any element of the algebra g^ formally has the same form as the associated 
element from g^, but in fact is a vector field in the different space augmented with the arbitrary 
elements. 




elements, belongs to G^. 



□ 
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Proof. Consider arbitrary vector fields Qq € g*^ and Q G . Denote the trivial prolongation 
of Qo to the arbitrary elements by Qq, so Oo £ S^- It is necessary only to prove that (jo £ 0"" and 
[QjQo] S g'^. Let Go = {'7o(e) = exp(e(5o)} and G = {T(e) = exp(eQ)} be local one-parameter 
transformation groups associated with Qq and Q, respectively. As Gq is a subgroup of the 
vector field Qq belongs to g^. 

For each sufficiently small e define the composition T{e) = TQ{—\/e)T{—\/e)TQ{\/e)T{\/e) 
and consider the vector field 



Tie) 

e=0+ 



which coincides with [Q,Qo], see e.g. [25, Theorem 1.33]. As both T{—\/e)To{\/e)T{\/e) and 
7o(— belong to Gq (cf. Proposition 3), the transformation T'(e) also is an element of Gq. 
Therefore, Q G g". □ 

As the kernel is included in the maximal Lie invariance algebra of any equation from the class, 
we should classify only subalgebras of the equivalence algebra that contain the ideal associated 
with the kernel. 

Example 1. In general, the kernel g^ is not necessarily an ideal of the maximal Lie invariance 
algebra Qg for each 9 G S. Indeed, consider the class of (1 + l)-dimensional nonlinear diffusion 
equations of the general form ut = {F{u)ux)x, where F ^ 0, cf. the introduction. The kernel of 
this class and the maximal Lie invariance algebra of the diffusion equation with F = u~^^^ are 
gn = (^Q^^ Q^^ 2tdt + xdx) and gi = {dt, d^, 2tdt + xd^, Atdt + 3udu, x'^d^ - 3xudu) [1, 27, 28], 
respectively. At the same time, the kernel q'^ is not an ideal of gi, [g^,gi] ^ g*^, since 

[dx,x'^dx — Sxudy] = 2xdx — Sxudu g^. 

Note that the class of diffusion equations is semi-normalized (see [30, 32] for the definition of 
semi-normalization) but not normalized in the usual sense. 

Corollary 2. If the class C\s is normalized in the usual sense, the kernel algebra g'"' is an ideal 
of the maximal Lie invariance algebra Qg for each 9 £ S. 

Proof. We fix an arbitrary element 9^ S. Denote by Qgo the maximal subalgebra of g"^ such that 
the algebraic equation 9 = 9'^{x, u^r)) is invariant with respect to it. This subalgebra necessarily 
contains the trivial prolongation g*^ of the kernel algebra q'^ to the arbitrary elements. Thus, we 
have that g^ C g^o C g~ and, in view of Corollary 1, g^ is an ideal in g~. Therefore, g^ is an 
ideal in Qgo. As the class C\s is normalized in the usual sense, the projection of g^o to the space 
of the variables (x, u) coincides with the maximal Lie invariance algebra g^o of the equation Cgo. 
By the construction, the projection of g^ to the space of the variables {x,u) coincides with g^. 
Hence g^ is an ideal in g^o . □ 

Often the equivalence algebra can be represented as a semi-direct sum of the ideal associated 
with the kernel and a certain subalgebra. To obtain preliminary group classification in this 
case, we in fact need to classify only inequivalent subalgebras of the complement of the kernel 
ideal. Projections of these subalgebras to the space of equation variables will give all possible 
inequivalent extensions of the kernel. 

Example 2. We present a class of differential equations for which the above representation is 
not possible. This is the class of (1 + l)-dimensional linear second order homogeneous evolution 
equations which has the general form 

Ut = A{t, x)uxx + B{t, x)ux + C(t, x)u, (2) 
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where A = A{t,x), B = B{t,x) and C = C{t,x) are arbitrary smooth functions, A ^ 0. The 
kernel Lie algebra of class (2) is = {udu)- Its equivalence algebra g~ is spanned by operators 
of the form 

rdt + Cdx + V^udu + 

(2^. - Tt)dA + ((e. - rt)B - 2r,lA - it)dB + {vl " Mx " Bvl - CCt)dc, 

where r = r(t), ^ = ^(t, x) and t]^ = r]^{t, x) are arbitrary smooth functions of their arguments. 
The kernel can be identified with the ideal g'"' of g~, generated by the vector field udu, which 
is assumed now to act in the space of variables and arbitrary elements. Moreover, this vector 
field commutes with all elements of g~. At the same time we have [dt,tudu] = udu- Therefore 
the algebra g" cannot be represented as a semi-direct sum of g^ and a subalgebra. 

3 Determining equations of Lie symmetries 

The method of computing Lie symmetries is classical and can be found in all textbooks on this 
subject, see, e.g. [6, 25, 28]. Owing to its algorithmic nature, it was implemented in a number 
of symbolic computation programs [7, 8, 10, 35]. For an equation A = from the class (1), the 
condition of infinitesimal invariance with respect to a vector field 

Q = T{t, X, u)dt + i{t, X, u)dx + X, u)du 

has the form Q^^-'Aj^^Q = 0, i.e., 

g(2)A = rf - Cfxul - r]fuul - 2fuxV'' - ^QxUxx - WuUxx - QV'"'' = (3) 
wherever A = 0. Here Q^^^ is the second prolongation of the vector field Q, 

where the coefficients can be determined by using the general prolongation formula. In (3) we 
only need the coefficients r/*, rf and rf^. They read [25, 28] 

if = Bt{ri- TUt - ^Ux)+ TUtt + Cutx, 

T]'' = Dx{r]-TUt-(,Ux)+TUtx+CUxx, (5) 
^xx ^ -p2^^ _ _ ^^^-j _^ ^^^^^ iUxxx, 

where and Dx denote the operators of total differentiation with respect to t and x, respectively, 

Dt = dt + utdu + uttdut + utxdu^ H , D^. = 9^. + Uxdu + Utxdut + Uxxdu^ H • 

Upon plugging the coefficients (5) into the infinitesimal invariance condition (3), we obtain 
the following equation 

DtT] - utBtT - UxBtC - ifxul - i]fuul - 2fux(Dxr] - utB^T - UxD^O - COxUxx - 

2 2 2 y^) 

VduUxx - gi^xV - Ut^xT - ""xD^^ - 2utxT)xT - 2UxxGxi) = 0. 

In order to constrain this equation on the manifold of equations (1), we set Ut = fu^ + gUxx- 
Then, splitting (6) with respect to the various derivatives of u we obtain the following overde- 
termined system of determining equations of Lie symmetries: 



UxUxt ■ 


Tu = 0, 


Uxt- 


Tx = 0, 


'^x'^xx • 


e« = o. 




fin + r]u- 2^x) + gvuu + Ux + vfu 


Uxx • 


gin - 2Cx) + Cgx + = o. 


Ux- 


6 + 2fr]x + g{2r]xu - Cxx) = 0, 


1: 


m - grixx = 0. 



(7) 
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As usual for classes of differential equations, the determining equations split into a part not 
involving the arbitrary elements and a part explicitly involving them (the classifying part). 
In the present case, the first three equations do not involve / and g and can therefore be 
integrated immediately. They give r = r(t) and ^ = ^{t,x), i.e. the symmetry transformations 
are projectable and transformations of t only depend on t. 

The remaining four equations form the system of classifying equations. In the case of arbi- 
trariness of the functions / and g, we can further split system (7) with respect to derivatives of / 
and g. This yields the kernel of maximal Lie invariance algebras, which gives rise to those sym- 
metry transformations that are admitted for all elements of the class of equations (1). Splitting 
yields 

e = 7/ = 0, Tt = 0, 

i.e. the kernel algebra is generated solely by the operator dt, 0^ = {dt)- That is, for arbitrary 
values of / and g, the only symmetry admitted by equations of the class (1) is the time translation 
symmetry (t, x, u) {t + e,x,u), e £ M. 



4 The equivalence algebra 

In order to investigate inequivalent cases of symmetry extensions of the kernel algebra q'^, the 
equivalence algebra (group) must be computed. Unfortunately, the equivalence algebra presented 
in [20] is not correct. It can easily be checked that their operator h{x)du does not generate 
an equivalence transformation for general values of h. Similarly, also their operator a{x)dx + 
2fa{x)df + ga'{x)dg cannot generate equivalence transformations for arbitrary values of a. The 
problem indeed is that their infinitesimal invariance condition for equivalence transformations 
is incorrect. This is why it is necessary to re-derive the equivalence algebra for the class of 
equations (1) here. 

Theorem 1. The equivalence algebra 0~ of the class of equations (1) is generated by the fol- 
lowing operators: 

dt, dx, V' = tdt- fdf-gdg, = xdx + 2fdf + 2gdg, 

g{h) = hdu - {huf + huug)df, 
where h = h{u) is an arbitrary smooth function of u. 

Proof. The proof is done using infinitesimal methods. We seek for operators of the form 

Y = Tdt + + r^du + ifdf + edg 

that generate continuous equivalence transformations, where r, ^ and ij are functions of the 
variables t, x and u, whereas and 9 are regarded as functions of t, x, u, f and g. That is, 
we aim to determine the usual equivalence algebra rather than some generalized equivalence 
algebra [19, 32]. The class of equations (1) must be augmented with the auxiliary system 

Si := ft = 0, S2 := gt = 0. (9) 

The complete auxiliary system should also include the conditions that the arbitrary elements / 
and g do not depend on nonzero order derivatives of u. However, these conditions already are 
implicitly taken into account by the supposition that the coefficients of Y does not involve these 
derivatives. 

The joint invariance condition then reads 

yA|_^=0, YSi\j^=0, YS2\j^=0, (10) 
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where M denotes the joint system of the equations A = 0, 5i = and 5*2 = 0, 

and Q^^) is defined by (4). The coefficients and 0* can be obtained by the first prolongation 
considering (t, x, u) and (/, g) as independent and dependent variables, respectively, 

=t)t{ip- Tft - Cfx - rjfu) + rftt + Cftx + r]ftu, 
6** = t)t{ip- TQt - CQx -r]gu) + TQu + igtx + i]9tu, 

where Dt = 9j + ftdj + gtdg + • • • is the corresponding operator of total differentiation with 
respect to t. In view of the auxiliary system (9), the total derivative operator reduces to the 
partial derivative, i.e. Dt = dt- 

The second and the third conditions from (10) then imply that 

n - Ctfx - vtfu = 0, Ot- Ctgx - vt9u = 0. 

Since these equations should be satisfied for all values of the arbitrary elements / and g, we can 
split with respect to the derivatives fx, fu, dx and gu to obtain that 

Vt = dt = Ct = Vt = 0. 
It remains to investigate the first condition in (10). In detail, it reads 

?7* - 2fuxV'' - ful - gr/^'^' - Ouxx = 0, 
or, after expanding, 

DtT] - utBtT - UxBtS, - 2fuxiBxT] - utBxT - UxDxO - - 
giBlr] - utBlr - Ux^l^ - 2utx^xT - 2uxx^xO - Ouxx = 0. 

We now split this equation with respect to the derivatives of u similar as done in the course of 
deriving the determining equations of Lie symmetries. The splitting with respect to utx implies 
that r = T{t). Splitting with respect to UxUxx, we derive that ^ = ^(x). These conditions already 
simplifies the above invariance condition substantially. Collecting coefficients of the remaining 
monomials of derivatives leads to 

Uxx- = {2^x - Tt)g, 

ul- ^ = {2ix -Tt- r]u)f - rjuug, 

Ux ■ 2r]xf + 2r]xug - ixxQ = 0, 

1 : 9'nxx = 0. 

In view (ft = Of = £,t = Vt = 0, the general solution of this system is 

T = cit + C2, ^ = C3X + C4, ri = h{u), 

99 = (2c3 - ci - hu)f - KuQ, = (2c3 - ci)g, 

where ci, . . . , C4 are arbitrary constants and h is an arbitrary smooth function of u. 

This completes the proof of the theorem. □ 

The equivalence algebra 0~ can be represented in several ways, which are important for differ- 
ent purposes. The representation crucial for the present case is that = {dt)^{'D^, P*, dx, G{h)), 
see Remark 3 for further details. Another natural representation is 0~ = {dt € 2?*) © {dx S T>^) © 
{Q{h)). This representation implies that g"" is the direct sum of a finite-dimensional and an 
infinite-dimensional parts. This representation is helpful for the determination of the adjoint 
actions, see Section 7. 
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5 The equivalence group 



In the previous section we have determined the equivalence algebra of the class (1) using in- 
finitesimal techniques. In order to obtain the complete point equivalence group (including also 
discrete transformations), the direct method should be applied. For the sake of completeness 
and illustration we present the corresponding computations here. 

Theorem 2. The equivalence group of the class of equations (1) is formed by the transfor- 
mations 

i=Ait + Ao, x = Bix + Bq, u = U{u), 
~ Bl ( Uuu \ ^ Bf 

where Aq, Ai, Bq, Bi U is an arbitrary smooth function of u and AiBiUu ^ 0. 

Proof. We begin with a preliminary description of admissible transformations of the class (1). 
In other words, we derive determining equations for point transformations that map a fixed 
equation from the class (1) to an equation from the same class. As (1) defines a subclass of 
(1 + l)-dimensional evolution equations, we at once know that the transformation component 
of t depends only on t, see e.g. [13, 18]. Moreover, each equation from the class (1) belongs to 
the class of second-order quasi-linear evolution equations having the form ut = F{t, x, u)uxx + 
G{t, X, u, Ux)- Hence in view of Lemma 1 of [12] the transformation component of x depends only 
on t and x. That is, the transformations of variables will be of the form t = T{t), x = X(t,x), 
u = U{t,x,u), where TfXxUu ^ 0. The transformed derivatives then read 

Ui = ^ (BtU - ^^xU^ , Ux = -^^xU, Uxx = (^D^) ^■ 

Substituting these derivatives into the transformed form of (1) and taking into account the 
initial form (1), we obtain 



^(^t^t-^Dx.f/J+^(/<. + 5n..) = /(^^J +~g[j^J^xj U, (11) 

where / = f{X, U) and g = g{X, U). Splitting equation (11) with respect to Uxx and Ux yields 
X ^ 

Uxx'- g = ^^g, (12a) 



t 

ij ~ J J 

Equation (12a) defines the transformation rule for the arbitrary element g. Substituting equa- 
tion (12a) into (12b) leads to the transformation rule for the arbitrary element /, 

^ ~ TtUj TtUl 

In general, system (12) forms the determining equations for admissible transformations of the 
class (1). Unfortunately, this system is too difficult to be integrated since there are a lot of 



11 



different cases of its solution depending on specific values of the arbitrary elements. However, 
this system allows to easily determine the equivalence group. For this aim, we can split equa- 
tions (12c) and (12d) with respect to / and g. This gives at once Xt = X^x = Ux = Ut = since 
Uu ^ 0. Furthermore, differentiating the first equation of system (12) with respect to t leads to 
the final restriction Tu = 0. Solving these determining equations for equivalence transformations 
completes the proof of the theorem. □ 

Corollary 3. A complete set of discrete equivalence transformations in the group , which 
are independent up to their composition and composition with continuous transformations are 
exhausted by the three transformations of alternating signs 

It- {t,x,u,f,g)^{-t,x,u,-f,-g), 
Ix- {t,x,u,f,g)^{t,-x,u,f,g), 
lu- {t,x,uj,g) ^ {t,x,-u,-f,g). 

The equation (1) with the specific value 6*0 = if,g) = (— 4/3n~^/^, ti~^/'^) admits the Lie 
symmetry operator x'^dx — Sxudu- The transformations from the corresponding one-parameter 
transformation group belong to T{6q,6q). As the associated admissible transformations are not 
induced by elements of the equivalence group G~ of the class (1), this class is not normalized. 
Similar assertions are true for the potential Burgers equations {f = g = 1), linear equations 
from the class (1) (/ = 0, = 0), etc. As system (12) is too complicated and the equivalence 
group G~ is quite narrow in comparison with the class (1) (the transformations from G~ are 
parameterized by four constants and only a singe function of one argument and, at the same 
time, the tuple of arbitrary elements consists of two functions of two arguments), this justifies 
why preliminary group classification is well suited for the class of equations (1). 

Remark 2. It is not possible to simplify the general equation from the class (1) by equivalence 
transformations. The interesting particular case of simplification by equivalence transformations 
is given by equations of the form (1) with / proportional to g. If / = eg, where c is a nonzero 
constant, then the corresponding equation of the form (1) is mapped by the transformation 

i=t, i = x, u = e'''' (13) 

to the equation of the same form with / = and g = g{x, Inn). 



6 Classification of subalgebras 

In order to carry out preliminary group classification, it is necessary to derive an optimal list 
of inequivalent subalgebras. In the existing literature on the subject, usually only subalgebras 
of a certain finite-dimensional subalgebra of the equivalence algebra are classified up to inner 
automorphisms of this subalgebra. This restriction is, however, not necessary in the present 
case, although this is done in [20]. Furthermore, it should be noted that in [20] an algebra was 
chosen for preliminary group classification, which is not related to the corresponding equivalence 
algebra that was derived. 

To classify subalgebras of a Lie algebra of vector fields, it is necessary to know the adjoint 
action of the corresponding transformation (pseudo) group on this algebra. There exist two dif- 
ferent methods for the computation of the adjoint action. The first method employs information 
on the structure of the Lie algebra and is more suitable in the finite-dimensional case although 
it also works for certain infinite-dimensional algebras [5, ?, 9, 31]. The adjoint action of a 
one-parameter Lie group generated by an element v of the Lie algebra on this algebra can be 
determined either from the Lie series 

oo ^ 

w(e) = Ad(e^-)wo := V ^{v", wq}, 

n=0 
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where {v'^,wo} := wo,{v",wo} := (— l)"'[v, {v*^ ^,wo}], or, by solving the Cauchy problem 
^ = [w,v], w(0) = Wo, 

see [25] for more details. 

Following the second method, it is necessary only to calculate the actions of push- forwards 
of the transformations from the (pseudo) group on generating vector fields of the algebra. While 
the first method involves only the abstract structure of Lie algebras and therefore at once gives 
results for the whole class of isomorphic algebras, the second method relies on the specific 
realization of the Lie algebra by vector fields. At the same time, the second method works more 
properly in the infinite-dimensional case. 

We now derive the optimal lists of one- and two-dimensional subalgebras for the entire equiv- 
alence algebra 0~. 

The nonzero commutation relations of generating elements (8) of g~ are 

[d^,v^] = d„ [dt,v'] = dt, [g{h^),gih^)] = G{h^hl - h^hi). 

The nonidentical adjoint actions related to generating elements of 0~ and computed using the 
first method are 

Ad(e"^-)P^ = - ed^, M{e^'^^)d^ = e^d^, Ad(e^^('*'))g(/i2) = g{h^), 
Ad(e=^*)^* = - edt, Ad{e'^')dt = e'dt, 

where h?'{u,e) = h'^{H^{u, —e))/H^{u, —e) and {u = H^{u,e)} is the one-parameter transfor- 
mation group generated by the projection of the operator g{h^) to the space of the variable u, 
i.e. = h}{H^) and H^{u,0) = u. Although the four adjoint actions related to the finite- 
dimensional part of g'^ are suitable to be applied to the classification, there arises an inconve- 
nience with the adjoint action Ad(e^^'^^ ^) owing to problems with proving the existence of the 
required function h} . 

This is why, in what follows we use the adjoint action of the entire equivalence group C on 
the equivalence algebra , calculated by the second method. Any transformation T from G~ 
can be represented, for convenience, as a composition 

T = 'l'{Ao)7%Bo)V\Ai)D'-{Bi)^{U), 

cf. Theorem 2, where 



(14) 

gUuu/ul 

are translations with respect to t and x, scalings with respect to t and x and an arbitrary trans- 
formation of n, respectively, and AiBiUu ^ 0. Transformations of each of the above kinds form 
a subgroup of The last three subgroups contain the discrete transformations It, Ix and /„, 
respectively. Namely, It = 1), Ix = 2)^(— 1) and = ^{—u). As a result, additionally to 
avoiding the above problems with the existence of required values of functional parameters, in 
this way we at once include discrete equivalence transformations in the classification procedure. 

The nonidentical actions of push- forwards of transformations (14) on generating elements of 
0~ are exhausted by the followings: 

7-{Bo)V- = V- - Bodx, D-{Bi)dx = BA, %{U)G{h) = g{h{U)/Uu), 
7i{Ao)V' = - Aodt, 'Di{Ai)dt = Aidt, 

where the function U = U{u) is the inverse of U. 



T*(^o): 


i = 


t + Ao, 


X 


= X, 


u = 


u. 


9 


= 9, 


f = 


/, 


T-(So): 


i = 


t, 


X 


= x + Bo, 


u = 


u. 


9 


= 9, 


f = 


/, 


V\Ai): 


i = 


Alt, 


X 


= X, 


u = 


u. 


9 


= A^'g, 


f = 


Ai'f, 


2)^(5i): 


i = 


t, 


X 


= Bix, 


u = 


u. 


9 


= Bh, 


f = 


Bff, 


9{U): 


i = 


t, 


X 


= X, 


u = 


U{u), 


9 


= 9, 


f = 


f/Uu- 
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Remark 3. The kernel algebra generated by df is an ideal in the equivalence algebra which 
has the structure £|~ = (dt) € {'D^,'D^,dx,G{h)). Hence the classification of subalgebras of 
can be reduced to the classification of subalgebras of the algebra g~g = (P^, P*, dx,G{h)), which 
is the "essential" part of This will yield the possible Lie invariance algebra extensions of 
the kernel algebra obtainable by preliminary group classification. Moreover, the push-forwards 
of translations and scalings with respect to t should not be applied under the classification of 
subalgebras. 

Theorem 3. An optimal list of one- dimensional subalgebras of the algebra g^g is exhausted by 
the algebras 

^V- + aV' -0(8)), {v'+~5dx-G{5)), {d,-g{5)), {0(1)), (15) 
where a G M and 5,5 ^ {0, 1}. 

Proof. We use the approach for the classification of subalgebras that is outlined in [25] . We start 
with the most general form of an element of the algebra g~g, 

vi = a|P^' + + aldx + G{h^), 

where the constants a\, a\, a\ and the function = h^{u) are arbitrary but fixed, and sim- 
plify it as much as possible by means of push-forwards of transformations from the equivalence 
group In the case /i^ 7^ the function-parameter h^ can be set to —1 by usage of S*(f^) 
with the inverse [/ to a solution U = U{u) of the equation [/„ = —h^{U). In other words, up to 
G~-equivalence we can always assume that -h^ = 5 £ {0,1}. 

If a\ 7^ 0, we scale to set a{ = 1 and use the push- forward of a 7^{Bq) to set = 0. The 
notation a = a\ leads to the first subalgebra in the list (15). 

If a} = and a\ ^ 0, we set 02 = 1 by scaling and use D%{Bi) with certain Bi to set 
a\ = —6, where 5 G {0, 1}. This gives the second listed subalgebra. 

In the remaining case a\ = = we obtain the two last subalgebras from the list (15) under 
the assumptions Og 7^ and ag = 0, respectively, since the nonvanishing value of is set to be 
equal to 1 by scaling and the condition ag = necessarily implies that ^ and hence, up 
to G~-equivalence, = 1. □ 

Theorem 4. An optimal list of two-dimensional subalgebras of the algebra g^g reads 
ip'^ -g{5),V' -g{5)), (P" + aP* + g(n),5, -g(l)), (p- + aP* - ^(5), 9,) 

{V' -g{5),dx-Q{5)), + av' + bg{u),g{i)), {v' - 5dx + bg{u),g{i)), (i6) 
{dx-5g{u),g{i)), {g{i),g{u)), 

where a, b, 5, 5 and 5 are constants, and we can assume that 5,5 £ {0,1}, 5 £ if 5 = 1 and 
5 G {0,1} if 5 = 0. 

Proof. The proof of the above theorem is similar to those in the one-dimensional case, see a de- 
tailed explanation and other examples in [4, Chapter 7]. We start with two linearly independent 
copies of the most general element of g~, 

= alV^ + alv' + aldx + g{h^), 
v2 = afP^ + a^V' + a% + g{h^), 

and simplify them as much as possible by means of adjoint actions and nondegenerate linear 
combining. The additional complication concerns taking into account that the elements 
and should form a basis of a Lie algebra, i.e., their commutator should lie in their span, 
[v^, v^] G (v^, v^). Usually this places further restrictions on the admitted form of the elements. 
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To simply describe the conditions defining the different cases of the classification of two- 
dimensional subalgebras, we introduce the matrix notation 



( 



a 



a 



1 

'Ml 
,2 
'Ml 



a' 



a 



,1 

'Mn 
,2 

'Mn 



) 



where fii G {1,2,3} and n ^ 3. In what follows, the right hand side of a matrix equation 
^Mi - Mn = or a matrix inequality A^-^...^^ 7^ is the zero matrix of the appropriate dimension. 
In the course of classification, we should investigate two principal cases. 

1. rank(j4i23) = 2. This is the first cases which is partitioned into the three subcases 

(a) det^i2 / 0; (b) det A12 = 0, det A13 / 0; (c) det^i2 = 0, detvlis = 0. 

In the last subcase we necessarily have det A23 7^ 0. By means of a change of the basis we at 
first set A12 = E, A13 = E and A23 = E, respectively. Here E is the 2x2 identity matrix. 
If the new /i^ is nonvanishing, we set = -1 using S*(f^) with the inverse C/ to a solution 
of the equation C/„ = —h?{lJ). In other words, up to G~ -equivalence we can always assume 
that = 5 G {0,1}. We also set G {-1,0} in a similar way if h?' = 0. Specifically, in 
subcase (a) we further use the push- forward of T'^{a\) to set a\ = 0. As the resulting operators 
should commute, we derive that a| = and h\ = 0. This case hence leads to the first subalgebra 
from the list (16). In subcase (b) we re-denote a\ by a. Under the assumptions made, the 
commutator [v^,v^] equals — V2. Therefore, the condition /i2 = —1 implies that h\ = 1, i.e. we 
can set = u using a change of the basis and the push- forward of 7^{Bq) with certain Bq. This 
gives the second subalgebra from the list (16). If /i2 = 0, we obtain the third subalgebra. In 
subcase (c), the corresponding subalgebra is commutative and hence /i^ = 0. Applying a scaling 
of V2 and the push- forward of D^{Bi) with certain Bi, we simultaneously set h^^h? G {—1,0} 
and hence construct the fourth listed subalgebra. 

2. rank(Ai23) ^ 1. Up to a change of the basis, we can assume that of = a| = a| = and 
hence ^ 0, i.e., analogously to the previous case we can set /i^ = 1 by some S*(C^)- Then up 
to a linear combining of and the commutation condition [v^, v^] G (v^,v^) implies that 
h\ = b = const and, therefore, we can set = bu. The four last algebras from the list (16) 
represent the subcases 

(a) a\ / 0; (b) a\ = 0,4/ 0; (c) a\ = 0, ^ = 0, a\ + 0; (d) a\ = a\ = a\ = 0, 

in which by a scaling of vi we can set a| = 1, 02 = 1, 03 = 1 and 6 = 1, respectively. For the 
basis elements to have the appropriate canonical form, we should additionally set = by 
some TJ(i3o) and re-denote by a in subcase (a) and also set 03 G {—1, 0} by some 2)J(-Bi) in 
subcase (b) and b G {—1,0} by some 1)%{B\) and a scaling of vi in subcase (c). 

This completes the proof of the theorem. □ 

Note that all except the last subalgebras from the lists (15) and (16) represent parameterized 
classes of subalgebras rather than single subalgebras. 

7 Preliminary group classification 

Based on Proposition 1 and the above classification of subalgebras, we can obtain the extensions 
of the kernel algebra (9f) within the class (1) by projections of inequivalent one- and two- 
dimensional subalgebras of the equivalence algebra g~ to the space of variables {t,x,u). As 
a first step, for each of the subalgebras we solve the associated invariant surface condition 
for (/, 5), namely, the system of equations ^fx + rjfu = ^p, S,9x + V9u = ^1 where the operator 
rdt + S,9x + V^u + V'df + Odg runs through a basis of the subalgebra. 
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In Tables 1 and 2 we collect the general solutions of the invariant surface condition for (/, g) 
(or, in other words, the entire subclass of the corresponding invariant equations), which is asso- 
ciated with the one- and two-dimensional subalgebras of 0~ listed in (15) and (16), respectively. 
In these tables, / and g are arbitrary functions of single arguments and ci and C2 are arbitrary 
constants such that g ^ and C2 ^ 0. 

Table 1: One-dimensional Lie symmetry extensions for class (1) related to q^. 



N 


/ 


9 


Additional operator 


1 


f{u + 6ln\x\)x^-'' 


g{u + 51n \x\)x^~"- 


atdt + xdx — 


2a 


f{u + 5x)e~^ 


g{u + 6x)e~^ 


tdt + dx- Sdu 


2b 


f{x)e^ 


5(x)e" 


tdt - du 


3 


f{u + 6x) 


g{u + 5x) 


dx - Sdu 


4 


m 







Table 2: Two-dimensional Lie symmetry extensions for class (1) related to g"". 



N 


/ 


9 


Additional operators 


1 




c^e^x^^^ 


xdx -5du,tdt- du 


2 


ci\u + x\^~^ 




atdt + xdx + udu, dx - du 


3a 






atdt + xdx - du, dx 


3b 




5(w) 


2tdt + xdx, dx 


4 






tdt - du, dx - 6du 


5 


Ci|x|2-«-^ 


C2 Xp~" 


atdt + xdx + budu, du 


6 




C2e^ 


tdt - dx + budu, du 


7 




C2 


dx - dudu, du 


8 





~9{x) 


du, udu 



The second algebra from the list of one-dimensional subalgebras (15) is associated with a 
symmetry extension of an equation from the class (1) if and only if at least one of its parameters 6 
and 6 does not vanish. In addition, to find the corresponding ansatzes for / and g it is necessary 
to consider different cases of values of the parameters. This is why this subalgebra leads to 
two cases (2a and 2b) of Table 1. Analogously, equations from the class (1) are invariant with 
respect to the projections of the first, fourth or sixth algebras from the list of two-dimensional 
subalgebras (16) if and only if 5 7^ 0, i.e., we can assume that 5 = 1. For the third algebra we 
should have either 6^0 (then we can again assume that 5 = 1) or (5, a) = (0, 2) that gives 
Cases 3a and 3b, respectively. 

There are several reasons why Tables 1 and 2 do not give a proper classification result. We 
present these reasons in the form of the following series of remarks. 

Remark 4. As the whole consideration is done up to G~-equivalence, we should additionally 
factorize the general solutions of the invariant surface conditions for / and g with respect to 
this equivalence. Using transformations from G~, in Table 2 we can set C2 = 1 (by scaling of t 
and alternating its sign) and, in Cases 5;,=0) ^b=o and 7^=0 1 ci = (by the transformation (13) 
with c = Ci/c2, cf. Remark 2). For the other values of the parameters b and 6 in these cases, 
the constant ci can be assumed, up to G~-equivalence, to belong to {0, 1}. If a 7^ 2 in Case 3, 
we can scale the value 2 — a to 1. 



16 



Remark 5. Extensions presented in Tables 1 and 2 are not necessarily maximal even for the 
general values of the parameter- functions / and g or the constant parameters ci and C2- It 
lies in the nature of preliminary group classification that equations can admit operators which 
are not projections of operators of the equivalence algebra. For example, in the last case of 
Table 2 any corresponding equation is linear and therefore admits an infinite-dimensional Lie 
invariance algebra including also the operators of the form ip(t,x)du, where if runs through 
the set of solutions of the equation under consideration. (Of course, for certain values of g 
this equation possesses an even wider Lie invariance algebra, cf. [15, 28].) A similar remark is 
true for Case 5f,=o> (resp. Case 6;,=0) resp. Case 7^=0) of Table 2 since each of the equations 
corresponding to this case is reduced by an equivalence transformation to the linear equation 
with / = and g = (resp. g = e^, resp. g = 1), cf. Remark 4. 

Remark 6. What is more essential is that presented extensions are not maximal even among 
extensions related to subalgebras of g~. In particular, Case 3^=0 of Table 1 coincides by the 
arbitrary elements with Case 3b of Table 2 and hence should be excluded from the extension 
list. Within Table 2, if a = 2 the arbitrary elements in Cases 3a and 5f,=o coincide with those 
of Case 7^=0. Hence in Case 7^=0 we have the additional operator 2tdt + xdx induced by 
the operator + 2D*. The algebra presented in Case 3aj^2 is also not maximal, cf. Case 1 of 
Table 3. Cases 4, 6 and 7 admit additional extensions by the operator udu if ci = or e~^^'^^^^du 
if ci 7^ and 6 = (resp. 5 = 0), owing to the connection of these cases with Case 8 via the 
transformation (13). 

Remark 7. An effect of the lack of maximality of extensions is that under the simplification 
of the form of arbitrary elements by equivalence transformations the corresponding invariance 
algebra may be replaced a similar one. Thus, under setting ci = in Cases 5b=o, Qb=o and 7^=0 
the basis element du is replaced by udu- 

In order to complete the preliminary group classification of the class (1), we should at first con- 
struct the exhaustive list of G~-inequivalent subalgebras of whose projections to the space of 
the variables (t, x, n) are Lie invariance algebras of equations from the class (1). For convenience 
such subalgebras will be called appropriate. Then we should study the problem whether these 
subalgebras are maximal among the subalgebras with the same property for a certain subclass of 
the class (1). The majority of one- and two-dimensional subalgebras of g"" are appropriate. This 
is why for subalgebra dimensions one and two it is not too important whether all or only appro- 
priate subalgebras are classified but this is not the case for greater dimensions. As the arbitrary 
elements / and g depend on two arguments, the condition that the associated projection is a Lie 
invariance algebra of an equation from the class (1) is a strong restriction for subalgebras of g~ of 
dimension greater than two and even leads to the boundedness of dimension of such subalgebras. 

Let = {V\ Q{h)), where h runs through the set of smooth functions of u. For a subalgebra s 
of g~, we denote dims H g^" by m^. 

Lemma 1. D* s and ^ 2 for any appropriate subalgebra s o/g~. 

Proof. Let s be an appropriate subalgebra of g"". Then the system of invariant surface conditions 
associated with elements of s should have a solution {f^,g^) with g^ ^ 0. The invariant surface 
condition for g associated with the operator D* is 5 = that contradicts the auxiliary inequality 
£/ / 0. Hence V* s. 

In what follows, the indices i and j run from 1 to 3. Suppose that the subalgebra s contains 
at least three linearly independent elements from g^, v* = ^(/i*) + a'D*. The corresponding 
invariant surface conditions for g form the system h'^g^+a^g = 0. We consider it as a homogenous 
system of linear algebraic equations with respect to {gu,g)- This system should have a nonzero 
solution since g ^ 0. Therefore h^a^ — h^a^ = 0. In view of the linear independence of v^, 
and v^, this implies that all a* = and thus gu = 0. Now we interpret the system of invariant 
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surface conditions /i*/n + ^uf + Kiu9 ~ ^ / as a homogenous system of linear algebraic 
equations with respect to (/„, f,g). A.s g ^ 0, this system should possess a nonzero solution and 
hence the determinant of its matrix vanishes. At the same time, the determinant coincides with 
the Wronskian of the linearly independent functions h^, h? and /i^, which is not equal to zero. 
The contradiction obtained implies that ^2. □ 

Corollary 4. Any appropriate subalgebra 5 of g~ is of dimension not greater than four. 

Proof. The projection of any element from s \ gj" to (J)^,dx) should be nonzero. Therefore, 
dims ^ dim('D^, 5^:) + = 4. □ 

Corollary 5. s n gj" = s n {G{h)) for any appropriate subalgebra 5 of with = 2, where h 
runs through the set of smooth functions of u. 

Proof. As nis = 2, the subalgebra s contains two linearly independent elements from g^, v* = 
Q{h^) + 0*2?*, i = 1,2. Analogously to the proof of Lemma 1, we consider the system of the 
invariant surface conditions h^gu + a^g = for g associated with v* as a homogenous system 
of linear algebraic equations with respect to {gu,g), which has a nonzero solution since g ^ 0. 
Therefore, the determinant of its matrix equal zero, /i^a^ — h'^a^ = 0. In view of the linear 
independence of and v^, this implies that = = 0. □ 

As we have classified all one- and two-dimensional subalgebras of g~, it is enough to describe 
appropriate subalgebras only of dimensions greater than 2. 

Theorem 5. A complete list of -inequivalent appropriate subalgebras of is exhausted by 
the following subalgebras: 

{V^ + g(2), dx,V'- G{1)), (P-' + 2P* + 6g(n), a,., 

{v^ + av\g{i),g{u)), {dx-5v\g{i),g{u)), {v^ + 2v\dx,g{i),g{u)), 

where a, b and 5 are constants and we can assume that 6 G {0, 1}. 

Proof. Let s be an appropriate subalgebra of g~ and dims ^ 3. Then, > 0. Hence we should 
consider only the cases rUg = 1 and rug = 2. 

The condition = 1 means that the subalgebra s contains exactly one operator of the form 

= g{h^) + alV\ where / 0, in view of Lemma 1. By scaling of we can set = — 1 if 
a2 7^ 0. Moreover, as the function-parameter does not vanish it can be set to 1 upon using 
9*{U) with the inverse ^7 to a solution U = U{u) of the equation [/„ = h^{U). As a result we 
have two G~-inequivalent forms for v^: (i) = g{l) — f*, (ii) = g{\). The conditions 
dims ^ 3 and = 1 simultaneously imply that dims = 3. This is why we should have two 
more linearly independent operators of the form v* = a|P^ + a^I^* + a\dx + ^(/i*), i = 2, 3, from 
s \ g5" for which rank(a5^, a|)j=2,3 = 2, cf. the proof of Corollary 4. By linear combining of 
and we set a\ = a^ = 1 and a\ = a\ = 0. 

In subcase (i) we additionally subtract a^y^ from v* to obtain = in the new operator v*, 
i = 2, 3. The simplified form of and is = + g{h^) and = + g{h^), respectively. 
As s is a Lie algebra and {v^, v^, v'^} is a basis of s, any commutator of v's should lie in their 
linear span. This in particular implies that the operators and should commute with v^, 
which is equivalent to the conditions = and /if^ = 0. Then, the commutator [v^, v^] equals 
dx, which should belong to s. Therefore, = 0. The complete system of invariant surface 
conditions associated with {v^, v^, v^} has a solution with nonvanishing g if and only if h'^ = 2. 
As a result we obtain the first listed subalgebra. 

Analogously, in subcase (ii) we have [v^,v*] = ^(/i^) = where 6* = const, i = 2,3, 

i.e., up to linear combining of v* with v^, /i* = The condition [v^,v^] = dx ^ s yields 
that = fe^ = 0. In order to provide the requested compatibility of the entire system of the 
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associated invariant surface conditions with the inequality 5 7^ 0, we necessarily have = 2. 
Re-denoting 6^ = b, we recover the second subalgebra from the above list. 

If = 2, the subalgebra s contains two linearly independent operators v* = G{h^), i = 1,2. 
Similarly to case 2d of the proof of Theorem 4, we can assume up to G~ -equivalence and a change 
of the basis in (v"^, v^) that = 1 and /i^ = u. Consider any v = aiV^ + a2'D^ + a^dx + G{h) 
from the complement to {Q(l),Q{u)) in s. Lemma 1 implies that (01,03) 7^ (0,0). Therefore, 
[v*,v] G (v^,v^). The last condition is equivalent to hu,uhu—h G {l,u). Consequently, we 
obtain h € {l,u). Hence, up to linear combining with elements from {Q{l),0{u)), we can always 
assume that h = 0. In other words, the subalgebra s can be represented as a direct sum of the 
algebra {Q{l),G{u)) and a subalgebra of = {Vx,T>t,dx) ■ G~-inequivalent subalgebras of 
that do not contain the operator are exhausted by the algebras {T>^ + clD^), {dx — and 
{V^ + a'D^, dx), cf. the proofs of Theorems 3 and 4. In the last subalgebra, owing to the required 
compatibility of the system of invariant surface conditions associated with s we have a = 2. 

This completes the proof of the theorem. □ 

The symmetry extensions induced by subalgebras from Theorem 5 are collected in Table 3, 
where a is an arbitrary constant, a 7^ 2. Note that the extension induced by the second subalge- 
bra is not maximal among extensions related to g~. This is why we do not include it into Table 3. 
The general solution of the associated system of invariant surface conditions is / = ci and g = C2, 
where ci and C2 are arbitrary constants, C2 7^ 0. Such values of arbitrary elements correspond to 
the potential Burgers equation or the linear heat equation if ci 7^ or ci = 0, respectively. The 
linear heat equation is given by Case 4 of Table 3 and the potential Burgers equation, which 
additionally possesses the Lie symmetry operator e~^^^^^'^du induced by Q{e~^'^'^^^'^), is reduced 
to the same case by a transformation similar to (13), cf. Remark 6. Analogously, we should 
choose (5 = 1 in the fourth subalgebra for the associated extension to be maximal. 

Table 3: Lie symmetry extensions for class (1) related to g"^ of dimension greater than two. 



N 


/ 


9 


Additional operators 


1 






xdx + 2du, dx, tdt - du 


2 







atdt + xdx, du, udu 


3 







tdt - dx, du, udu 


4 





1 


2tdt + xdx, dx, du, udu 



Summing up the whole consideration of the present paper, we prove the following theorem: 

Theorem 6. The complete preliminary group classification of class (1) is split into Tables 1-3, 
where 5 ^ in Case 3 of Table 1 and in Table 2 we should globally set C2 = 1, exclude Case 3a 
and assume that 5 7^ —2 in Case 1, 5 ^ in Case 4, b ^ in Cases 5 and 6, and 6 = 1 in Case 7. 

Remark 8. Table 3 from [20], summing up the partial preliminary group classification of the 
class (1) therein, is incorrect. Neither are all of the equations listed really invariant under the 
operators presented in the table, nor are these operators proper additional operators in view of 
the kernel (dt). The main problem is that the basis element dt of the kernel is involved by linear 
combining to these additional operators which, moreover, are not linearly independent. 

The number of inequivalent cases to be investigated under the usage of the entire infinite- 
dimensional equivalence algebra g~ is rather small. This is due to the greater effectiveness of the 
adjoint action of the whole equivalence group, which allows for stronger simplifications under 
classification of inequivalent subalgebras. By using only a finite-dimensional subalgebra of g~ 
as usually done, the number of cases of extensions to be treated is generally greater. This is 
one more justification why it is favorable to use complete preliminary group classification rather 
than partial preliminary group classification. 
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8 Conclusion 



The main aim of this paper is a careful explanation of the technique of preliminary group clas- 
sification, its status in the picture of group classification, its benefits and its limitations. These 
points are those we consider to be mainly lacking so far. While preliminary group classification 
is generally attractive due to the relative simplicity of its algorithm, various of the results ob- 
tained by now using this approach have only little practical relevance, since they are presented 
without a detailed analysis of the class of differential equations. Moreover, in various papers 
only partial preliminary group classification was carried out, without indicating a sound physical 
justification for the chosen subalgebras of the respective equivalence algebras. Indeed, in some 
instances this choice might been motivated for the sake of pure mathematical convenience, which 
counteracts the initial aim of group classification of differential equations. 

In the present paper we substantially enhance the existing framework of preliminary group 
classification. We show that it is possible and convenient to treat subalgebras of the entire equiva- 
lence algebra even in the case if this algebra is infinite dimensional. This is the principal difference 
compared to existing works on the subject of preliminary group classification, in which the prob- 
lem is only partially solved by involving classification of subalgebras of a fixed finite-dimensional 
subalgebra of the equivalence algebra with respect to restricted adjoint actions. Furthermore, it 
is emphasized that only appropriate subalgebras satisfying certain properties should be classified. 

The algorithm of complete preliminary group classification can be summed up as follows: 

• Find the equivalence algebra 5~ and the equivalence group G~ of the class £1^ under 
consideration. 

• Classify appropriate subalgebras of g~ up to G~-equivalence, each of which satisfies the 
properties below: 

— It contains the kernel algebra oi C\s. 

— The associated system of invariant surface conditions with respect to the arbitrary 
elements is compatible. 

— It is the maximal subalgebra among all subalgebras of that have the same set of 
solutions for the associated systems of invariant surface conditions. 

• For each of the listed subalgebras, find the general solution of the associated system of 
invariant surface conditions with respect to the arbitrary elements. 

• Simplify these solutions using transformations from whose push-forwards to vector 
fields preserve the corresponding subalgebras of g~, i.e., these transformations lie in the 
normalizers of the corresponding subgroups of C. 

The systematic approach of complete preliminary group classification is exemplified with the 
class of generalized diffusion equation (1) that was recently attempted to be investigated in [20] 
using symmetry tools. Owing to the number of inconveniences of [20], we regard the class (1) 
as well-suited to explain the methodology of preliminary group classification. We use both the 
framework of the infinitesimal and the direct methods to derive the equivalence algebra and the 
equivalence group of the class (1). In addition, the direct method also allows us to obtain the 
classifying equations of admissible transformations. Similarly to the determining equations of 
Lie symmetries, these classifying equations are too difficult to be solved directly, which at once 
limits the chance to obtain a complete group classification of the class (1). 

It is important to indicate once more that the extensions of the kernel algebra constructed in 
this paper by using preliminary group classification are not necessarily maximal. That is, there 
are various equations in the class (1) which have the maximal Lie invariance algebras wider 
than the associated subalgebras of the equivalence algebra. This observation is another way of 
proving that the class (1) is not normalized. 
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